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Abstract 

We describe a numerical procedure to compute the so-called isospectral torus of finite 
gap sets, that is, the set of Jacobi matrices whose essential spectrum is composed of hnitely 
many intervals. We also study numerically the convergence of specific Jacobi matrices to 
their isospectral limit. We then extend the analysis to the definition and computation of an 
“isospectral torus” for Cantor sets in the family of Iterated Function Systems. This analysis 
is developed with the ultimate goal of attacking numerically the conjecture that the Jacobi 
matrices of I.F.S. measures supported on Cantor sets are asymptotically almost-periodic. 
keywords: Potential Theory Orthogonal Polynomials Iterated Function Systems Almost pe¬ 
riodic Jacobi matrices. 
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Dedicated to Ed Saff on his seventieth 


1 Introduction 

Harmonic analysis have always been one of Ed’s interests, and while working on problems of 
theoretical sophistication, he never disdained the applicative side of the theory. Years ago we 
had the occasion to discuss the usage of Szego polynomials in the determination of harmonic 
frequencies embedded in a real (or complex) signal (HT] [1] . Ed has been a pioneer of the 

mathematical analysis of this approach, that is now being applied to challenging problems like 
gravitational waves detection by a long-time friend of Ed’s, Daniel Bessis |58) . 

In this paper I want to discuss a theoretical problem in which harmonic analysis plays a 
fundamental role, and where numerical experiments are important, for lack of a definitive theory. 
About twenty years ago [33] I conjectured that the entries of the Jacobi matrix associated with 
the devil’s staircase measure are almost periodic (all terms to be defined precisely later on), and 
I provided the first numerical evidence to that aim. The conjecture was rather natural, since 
another example of singular continuous measure supported on a Cantor set, that is, Julia set 
measures, possess almost periodic Jacobi matrices, provably so at least in a certain range of 
parameters [B (nil [H [III, and the so-called almost periodic flu [3011131311111111] had just 
raged. Moreover, in the inverse direction, almost periodic Jacobi matrices, like those constructed 
according to the Fibonacci rule have a spectral measure supported on Cantor sets with a clear 
hierarchical structure [03 ■ 

Although this conjecture had been reported repeatedly afterwards [33111113, orthogonal 
polynomials on I.F.S. Cantor sets remained confined to the study of quantum evolution of systems 
with singular spectra, see e.g. |48j in Ed’s 60th anniversary volume and references therein. This 

* Computations for this paper were performed at the INFN computer center at Pisa. 
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topic, and the conjecture, were recently brought to evidence by two papers that were motivated 
on the one hand |33] by the search for approximate identities with Dirichlet kernels, but went 
much further, and on the other hand [39] by the infinite generalization of results for the so-called 
finite gap sets, that is, set of intervals, that opened new perspectives in this problem. It is 
therefore timely to return to this problem. Although concerned with this dated question, the 
present paper is not simply review, but it introduces new ideas and techniques, as the occasion 
of Ed’s seventieth deserves. As they say, happy continuation, Ed! 


2 Jacobi matrices on IFS sets: definitions and summary of 
the paper 


Orthogonal polynomials, {pj(/r; s)}jgN, of a positive Borel measure /i supported on a compact 
subset E of the real axis are defined in a straightforward way by the relation f pj (/i; s)pmip', s)dp,{s) 
Sjm-> where Sjm is the Kronecker delta. The well known three-terms recurrence relation 

spj{p; s) = bj+ipj+i{p] s) + ajPj{p] s) + bjPj-i{p] s), (1) 


initialized by 6o 
matrix J(/i): 


0 and p-i{p]s) = 0, po{p;s) = 1, can be formally encoded in the Jacobi 


J{p) := 



bi 

ai 


^2 


( 2 ) 


For compact support E the moment problem is determined |5], and the matrix J(/i) is in one- 
to-one relation with the measure p. 

In this paper we will study the Jacobi matrix of a variety of measures supported on set of 
intervals, defined by Iterated Function Systems, reviewed in Sect. |3l This construction produces 
a sequence of finite unions of intervals, {A„}„gN, converging to a Cantor set when n tends to 
infinity. Each i?„ is a compact set, that carries a unique equilibrium measure, A numerical 
procedure to compute this measure is described in Sect. 0] following [33]. This procedure yields 
the set of fundamental harmonic frequencies that are the basis of further developments. In sect. [5] 
a novel numerical method for the computation of the isospectral torus of the sets A" is described. 
The infinite Jacobi matrices so computed are then analyzed in their harmonic components in sect. 
[HI by using a refined technique outlined in Appendix. Validity of this technique is independently 
verified in the successive sect. 0 where we describe convergence to the isospectral torus of the 
Jacobi matrix elements of two significant I.F.S. examples. It might be that this is the first time 
that this convergence has been explicitly exhibited in a numerical example. 

The final two sections touch upon two challenging problems. First, in sect. [S] we propose a 
definition of the isospectral torus of infinite gap sets-Cantor sets generated by I.F.S. that natu¬ 
rally leads to numerical computations. Finally, in the closing sect. [3] we propose a framework to 
assess the validity of the conjecture mentioned in the introduction, that is, we discuss whether 
the Jacobi matrix of an I.F.S. balanced measure—like the devil’s staircase measure—is asymp¬ 
totically almost periodic. We study to what extent its matrix elements can be approximated 
by the trigonometric polynomials obtained by the harmonic analysis developed in the preceding 
sections. 


3 Attractors of Iterated Function Systems 

Let us first construct finite gap sets, converging to a Cantor set. To do this, we use the simplest 
version of Iterated Function Systems (IFS) [S3] [MllIIl [IH] ■ These are collections of contractive 
maps (pi : R —R, i = I,..., M, for which there exists a unique set A, called the attractor of the 
IFS, that solves the equation 

A= U MA)--=HA). (3) 
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The attractor A can also be seen the limit (in the Hausdorff metric) of the sequence i?" := 
<i>"(£’°), where is any non-empty compact set: A = lim„_).oo We now choose a finite 

number of affine maps of the form: 

(l>3is) = Sj{s - 7 j) +-fj, j = (4) 

where Sj are real numbers between zero and one, called contraction ratios, and "fj are real 
constants, that geometrically correspond to the fixed points of the maps. Furthermore, we choose 
the set as the interval E^ = [miny^, maxy^], and we also impose that (f>i{E'^) fl (j>j{E^) = 0 , 
for i ^ j. This gives rise to what is called a fully disconnected IFS, whose attractor ^ is a Cantor 
set. It follows from standard theory that E^ is the convex hull of A and that E'^ is the union of 
M" disjoint, closed intervals, that we will denote as i?" = [ai,/3i]: 

M" 

E^ = <^^{E°) = \J[a,,f3,]. (5) 

i=l 

In accordance with the spectral analysis of periodic solids, we will call the intervals Ef' at 
r.h.s. of eq. ([5]) bands at generation (or level) n. It appears from eq. ([5]) that all bands Ef 
at level n can be obtained by applying the transformations (fj, j = 1,... ,M, to the bands at 
generation n — 1 . As a consequence, the length of these intervals, as well as their sum, tend to 
zero geometrically and the attractor is a Cantor set of null Lebesgue measure. Complementary 
to bands, in the convex hull of A", are the so-called gaps, for which we will use the notation 

M"-! 

G":=|jGr= U (A,«.+i). ( 6 ) 

i 2=1 

4 Equilibrium measure for a set of intervals 

Crucial in the following is the role of the equilibrium measure for A", for an electrostatic charge, 
with a logarithmic law of repulsion. We now recall briefly the main facts of this theory . 

Let a any Borel probability measure, supported on A”. The potential V{a-,z), generated by tr 
at the point z in C, is 

V{a;z):=-f log |z - s| dCT(s). (7) 

J 

The electrostatic energy £{a) of the distribution a is the integral of V{a- z): 

8{a) := I V{<7-,u) da{u) = — j f log \u — s\ da{s)da{u). ( 8 ) 

JE" j j e^ 

The equilibrium measure ve^ associated with any compact domain A" is the unique measure 
that minimizes the energy £{a), when this latter is not identically infinite. This being easily 
proven in the case of a finite union of intervals, Cap(A") := defines the capacity of the 

set A", and the Green’s function can be written as 

5 (A”; z) = z) - log(Cap(F;")) = -^(^ 5 .; z) + £{ve^). (9) 

From a theoretical point of view, the solution of the equilibrium problem for a finite union 
of N intervals i = 1,..., fV, is well known. Let’s follow e.g. [69]. Define the polynomial 

Y{z), 

N 

Y{z)=W{z-(Xi){z-Pi), ( 10 ) 

2 = 1 

and its square root, ^JY (z), as the one which takes real values for z real and large. Also, let the 
real number Qi belong to the open interval {Pi, cci+i), for i = 1,..., iV — 1 {i.e. to the gap G"). 
Define Z{z) as the monic polynomial of degree TV — 1 with roots at all C^fs: 

N-l 

Z{z)=l[{z-C,). (11) 

2=1 
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Theorem 1 There exists a unique set of values i = 1,..., 1} that solve the set of coupled, 

non-linear equations 



Z{s) 


ds = 0, i = 1,..., N — 1. 


( 12 ) 


In terms of these values, one can write the equilibrium measure of the set E'^: for simplicity of 
notation, we denote it by := ve^- 


dv^{s) 


1 ^ 


l(s) 




\Zis)\ 

VWWl 


ds. 


(13) 


It is apparent from the previous equation that the measure is absolutely continuous with 
respect to the Lebesgue measure on Moreover, since the compact sets of positive capacity 
tend to A when n grows, standard theory (see e.g. 0) assures that 

Theorem 2 As n tends to infinity, tends weakly to the equilibrium measure on the Cantor 
set A. 


In two previous works [491151] , we have shown how to solve numerically the set of equations 
(|12l) in a stable way, for sets of intervals of large cardinality. In [51] we have also quantified the 
rate of this convergence in terms of the orthogonal polynomials z). 

Of particular relevance are the integrals of over the intervals composing S", ^”([ 0 ;^,/3i]). 
When these integrals are rational numbers, Peherestorfer has proven that there exists a 
strict-T polynomial on £1" , that is, a polynomial with oscillation properties mimicking, and 
extending, those of the classical Chebyshev polynomials on a single interval. This is notably the 
case when the I.F.S. maps defining E'^ are the inverse branches of a polynomial dynamics with 
a real Julia set (see [T3|i9i dg da ED] and [51] for a numerical investigation). On the contrary, 
this is not the typical case of intervals generated by affine maps of the kind (|4]). This fact has 
important consequences in the following. In fact, the harmonic analysis that we will develop 
features the fundamental harmonic frequencies: 


1 l^(^)l ,, 


(14) 


where is any point in the i — th gap Typically, for IPS affine maps, the values ujf, 

i = 1,..., N — 1 are not rationally related. 


5 Set of Intervals — The Isospectral torus 


The isospectral torus of Jacobi matrices associated with measures supported on a finite set of 
intervals arise naturally in the study of the asymptotics of orthogonal polynomials supported on 
these sets [sniisiiTiEgisi El m di uning ■ Therefore it has been extensively studied and many 
equivalent definitions shed light on its nature El EH- We want to adopt here a constructive 
definition mm that reads as follows. Let Y{z) be defined as above, let Z{^;z) = ni=i^('2“?i)) 
where now are arbitrary values in the gaps and let V be a polynomial constructed in such a 
way that the function 


M{z) 


X{z) - vT(i) 
Ziz) 


(15) 


behaves as 1/z for large z and X{fi) = ai^yY{fi), with |tTi| = 1. Clearly, the n — 1 numbers fi 
and the signs at label points in an n — 1 dimensional torus. It can be shown that 


Theorem 3 M{z) is the Markov function (a.k.a. Stieltjes transform) of a positive 


measure 9, 


M{z) 



d9{s) 
z — s' 


(16) 
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where 9 is given by the sum of an absolutely continuous measure and of a linear combination of 
atomic measures supported at the zeros of the function Z: 


d9{s) 


N 




2=1 




N 

ds + ^ ■ 

i=l 


(17) 


Observe that the constants A and Bi depend on the parameters and at via eq. (fTCl) . The 
absolutely continuous part contains a generalized Chebyshev measure of the second kind. 

Definition 1 The isospectral torus of a finite set of intervals is the set of Jacobi matrices asso¬ 
ciated with the measures 9 in eqs. mmn^- 

This theorem prompts for a numerical computation of such Jacobi matrices, adapting the 
techniques of m- In essence, the theory run as follows: 

• The absolutely continuous component in p7l) is written as the (weak) limit of a sequence 
of atomic measures (Lemma 2 in [51)1. constructed using the properties of Chebyshev mea¬ 
sures, of first kind for the equilibrium measure, of second kind for the present application. 

• The (truncated) Jacobi matrix of a finite sum of atomic measures is computed adapting and 
enhancing existent procedures. These techniques allow for a uniform control of convergence 
to the desired matrix elements: we can estimate max{|aj —dj|, \bj —bj\,j = 0,..., J}, where 
the bar indicate the numerically computed matrix elements, and reduce this distance within 
acceptable limits, even for large values of J. 

While our techniques are fully general, we need to choose a specific example for the numerical 
computations of this paper. At difference with Ref. m, where the ternary Cantor set has been 
chosen, we focus here on a different I.F.S., that does not enjoy of the useful symmetry properties 
of the former, and represents a more significant instance of the general case. This is the following 


Example 1 Let M = 2, Si = 0.34, 62 = 0.52, 71 = — 1, 72 = 1. 

In figure [1] we draw the initial parts of the sequences {aj}jLi, {bjfffi for a point in the 
isospectral torus of the set generated by the IFS of example [TJ Our technique allows for 
large-scale computations: in figure[ 2 ]an estimate of the infinity distance between the numerically 
computed elements and the true ones, i.e. ej = max{| 6 ; —bi\,l < j} is plotted versus j. We 
observe a less than linear growth of the error that permits to obtain extremely accurate values 
even for matrix sizes in the order of hundreds of thousands. 


6 Harmonic analysis of the isospectral torus of finite gap 
sets 

In the previous section we have shown how to compute numerically the entries of Jacobi matrices 
in the isospectral torus of a finite number of intervals. It must be remarked that the same 
objects can also be expressed in terms of theta functions but this procedure seems to be 

numerically hard even for a small number of intervals. To the contrary, the approach adopted 
herein is capable of dealing satisfactorily with as many as 2 ^® intervals with moderate computer 
resources m- Before proceeding further with numerical techniques, we need an analytic result: 

Theorem 4 ( |54| I There exist real analytic functions 4>, T on the torus [0,1]^ ^ such that the 
recurrence coefficients aj,bj of Jacobi matrices in the isospectral torus can be expressed as 

aj = + f), bj = -b (j)), (18) 

where uj is the n — 1 dimensional vector with components ujf in eq. 0 , and (j) and ijj are shift 
vectors that depend on the particular matrix. 
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Figure 1: Diagonal (lower data) and outdiagonal (upper data) entries of a Jacobi matrix in the 
isospectral torus of (IFS of Example [T]) . 



Figure 2: Error in the outdiagonal elements of the numerically computed Jacobi matrix of fig. [T] 
versus j. The straight line has unit slope. 
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Figure 3: Jacobi Abel function in the case of a set composed of three intervals. 

This result is the basis of our numerical approach to the isospectral torus. For instance, Figure 
[3] displays the function 'I'(ai,a 2 ) in a case composed of three intervals, i.e. two fundamental 
frequencies, a case that permits a graphical rendering in three dimensions. How to parameterize 
conveniently this function is the topic of our further development. 

Let us from now on focus on the sequence of the outdiagonal matrix Jacobi matrix elements 
that is more significant mathematically than the sister All results will nonethe¬ 

less apply also to this latter. Let us expand the analytic function in Fourier series. Indicating 
by a the vector of arguments oi,..., aAr_i, a A^ — 1 dimensional real vector, and by A: a A^ — 1 
dimensional integer vector, we can write 

<l>(a) = Ck exp (27rzk ■ a). (19) 

ke/C 

The set of indices K. is to be chosen so to reflect the reality of the function dt. Using this 
representation in eq. ra, one trivially finds that 

bj = -f V') = E Ck exp (27rzk •'(/)) exp (2n7rik ■ w). (20) 

kGK 

However simple, the above equation reflects a significant, well known fact: 

Lemma 1 The outdiagonal sequence {bj}^i of a Jacobi matrix in the isospectral torus is almost 
periodic, with frequency module generated by the fundamental frequencies df. 

In equation (j20ll . one observes the frequencies Wk := 27rk-a; and the phases Uk := 27rk •■!/). Recall 
that the function dt does not depend on the specihc matrix in the isospectral torus. Therefore, 
we discover that 

Lemma 2 The almost-periodic spectrum of the sequence {&j(^, tT)}jgN for any choice of^, a is 
characterized by the same set of frequencies; the coefficient of the frequency tuk is Ck exp (z'i/jk) 
and therefore it has a constant modulus, irrespective of the choice of a. 

The coefficient C(w) of an almost periodic sequence is defined formally as the limit 

C(a;)= hm (21) 
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Figure 4: Amplitudes ICk] for example[Il 


We can therefore restate our approach to the isospectral torus associated with the set A" as 
follows: 

Corollary 1 The Jacobi-Abel function 4', represented as in eq. m can be determined via a 
Fourier analysis of the sequence for any single set of values 

To be rendered effective, this corollary needs a refined Fourier analysis of the sequence 
{bj}je'N- This can be achieved in three steps. 

• We compute the frequencies Wk- For this, we solve the set of equations (na) and we compute 
the fundamental frequencies ujf, following [4^1^ . 

• Using this information, we consider the numerical problem of computing the coefficients 
(^(wk) in a more effective way than by the definition (l2T|l that is theoretically significant 
but numerically ineffective. A technique to this aim is described in Appendix, when the 
set of indices JC is replaced by a finite set. 

• Finally, we organize the infinite set of indices as an increasing sequence of finite, nested 
sets. 

The first two steps above being described in Appendix and in the references, let us concentrate 
on the third. A possibility in this respect is to use a metric in that includes /C, so to define 

ICl = /C n Bl 

where Bl is the ball of radius L centered at the origin. Convergence of the Fourier sequence (IIID, 
when )C is replaced by /Cl and L tends to infinity is a classical problem of harmonic analysis, see 
e.g. [25]. In Figure |4| we use the metric induced by the norm ||k|| = |fci|, we set L = 8 and 

we consider A", with n = 2. In this case the support of the measure is made of four intervals: 
there are three fundamental frequencies tvf and k is a three-dimensional vector. At the location 
of this vector in space, we display a color-coded point, in logarithmic scale, according to the 
amplitude ICk]. We observe an exponential decrease of the amplitudes as ||k|| increases, that 
supports the relevance of the metric adopted. Yet, at an accurate inspection, this decrease is not 
uniform, but this observation leads to further discovery. 
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Figure 5: Amplitudes |(7k| for E^, examplelU versus ki, when k = (0,..., fci, 0,...), i = 1,..., 7. 


In fact, let us plot the values of |Ck| along lines from the origin directed as the coordinate 
axes. That is, fix all ki values to zero except one, and raise this latter from 1 to L. We obtain 
fig. ®, now for E^. The behavior observed for n = 3 is typical, and consists of an exponential 
decay of amplitudes with ki. 

The decays starts off from the value of the amplitude of the principal harmonic w”, that for 
simplicity we denote \C^\. This harmonic is associated with the gap G" in the set A". It has 
been observed long ago in the theory of electronic states in solids, that an harmonic component 
in the potential of a tight binding Hamiltonian opens a gap in the spectrum of this latter. By 
perturbation theory arguments m it was established that the size of the gap is proportional to 
the amplitude of the harmonic component. We observe here the reverse behavior: proportionality 
is approximatively well verified, but it is the size of the gap that induces the amplitude of the 
principal harmonic: see fig. IHl 

Remark that data points in figure |6] refer to four different orders n. Also notice that some 
points at different values of n appear to be almost coincident. These observations are of relevance 
when considering the limit of infinite n, a thing that we shall do later on. 

To conclude this section, we show in fig. [3 the harmonic amplitudes Gk versus Wk obtained by 
our technique for n = 3. Remark two fundamental differences with respect to fig. 2 in [44] (or the 
analogous fig. 3.1 in [39)1: a difference in technique—fast Fourier was there employed, limiting 
resolution, and a quite evident difference in symmetry. We are now in the presence of a typical 
spectrum (although still at a low value of n) of an IFS system. Also to test the reliability of the 
data so displayed, in the next section we want to tackle a difficult problem, that of convergence 
to the isospectral torus. 

7 Convergence to the isospectral torus 

In this section we will study the convergence of Jacobi matrices of different measures supported 
on A" to their limit point in the isospectral torus. We will consider the following situation: 

Hypothesis 1 Suppose that is an absolutely continuous measure with respect to the Lebesgue 
measure supported on A". Also suppose that p satisfies a Szego condition on A". 
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Figure 6: Principal harmonic amplitudes [Cf | versus |G"|, the size of the i-th gap in FI”. Values 
for n = 2,..., 5 are indicated by different symbols, example [TJ The fitting line is g{x) = Ax, 
with a suitable constant A. 



Figure 7: Leading harmonic amplitudes Ck versus Wk, for n = 3, ||k|| < 6. 
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We denote by {aj(/i„)}jgN, {^j(/^n)}jeN the sequences of Jacobi matrix elements of fin- In this 
respect, the following theorem is notable 

Theorem 5 ( j57j ) Under hypothesisUlthe recurrence coefficients aj{fj,n) and bj{p,n) are asymp¬ 
totically convergent to the almost periodic form given by eq. m, where the functions $, and 
the vector uj are derived from if", while the shift vectors (f and ip depend on p. 

A comment is in order. The above theorem has been proven in m in a much more general 
setting than the one presented here. On the other hand, a strengthening of this result is possible 
for special classes of measures. Aptekarev [S] considered absolutely continuous measures on if" 
with densities on each interval if" of the kind p{s) = w{s){s — ai)°‘'{bi — being aiffii > —1 
and w{s) locally analytical and non-negative. He proved that distance from the asymptotic 
limit provided by eq. (HU is asymptotically 0{l/j). Suetin [Bl], by requiring that p{s) = 
w{s)/ \/|T(s)| with a strictly positive, holomorphic w, proves that convergence is 0((5"), with a 
suitable constant J < 1. We will now study numerically both cases. 

To define measures satisfying Aptekarev or Suetin conditions, we make recourse once more to 
the IFS construction. Let us associate to each map (fj in eq. o a positive weight tTj , — 1 

and let us define the operator T acting on continuous functions on R given by T/ = Ylj ° 4>j- 
Let T* be the dual operator of T, acting in the space of regular Borel probability measures. 
Then, the following theorem is standard 

Theorem 6 Let pq any Borel probability measure on = [rniny^-,maxy^]. Let pn = 

(r*)"^o- Then /i„ tends weakly and in Kantorovich distance to a measure p supported on the 
attractor A, that is called an IFS balanced measure and satisfies the fixed point condition T*p = 

p. 

The devil’s staircase/Cantor measure considered in the Introduction is precisely of this kind. 
Observe that Theorem [6] applies to any initial measure po- We therefore make two choices. In 
case a) we use the normalized Lebesgue measure on E°. In case b) po is the (suitably rescaled) 
orthogonality measure of Chebyshev polynomials of second kind. It then easily follows that 

Lemma 3 Let /i„ be defined as above, with the IFS maps 0 and initial measure po as in cases 
a), b). Then, for any n, pn is an absolutely continuous measure that verifies hypothesis\^ In 
case a) and b) Aptekarev condition is satisfied, while Suetin condition holds only for case b). 

Let us now apply a numerical technique for the stable computation of the Jacobi matrices 

, that applies to both cases a) and b) just discussed. This technique has been introduced in 
ref. [46] and is radically different from the one described in m and in Sect. [Hof this work. On 
the force of the above theorems, these matrices converge to a point in the isospectral torus. This 
raises a challenging numerical problem. 

Problem 1 Given a measure p satisfying hypothesis\^ and its Jacobi matrix J^, find the matrix 
in the isospectral torus of if" to which converges. 

We propose the following solution of this problem. In view of the harmonic analysis presented 
in Sect. [Bland of Lemma [5] we can reduce the previous problem to the determination of the shifts 
V'k, for k € /Cl. It is then possible from these latter to reconstruct the precise isospectral matrix, 
since amplitudes and frequencies are universal in the isospectral torus and can be determined in 
advance. We denote by and Jg^ measure and matrix in the isospectral torus that correspond 

lO pn . 

Again, let us consider for simplicity the single sequence The associated set of 

phases '0k can be computed combining Fourier analysis with a suitable extrapolation procedure 
outlined in the second part of the Appendix. Let us comment here on the results that can be 
obtained in this way. 

We start by examining case a). Fig. [8[ exhibits convergence of bj{pn) to the corresponding 
matrix element bj{9n) of the measure in the isospectral torus of if", obtained as described 
above. Convergence is slow enough to permit us to appreciate visually how the two sequences 
are getting close. 
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Figure 8: Outdiagonal entries of the Jacobi matrix of (circles) and of its limit 6 n in the 
isospectral torus (diamonds). Here, n = 4 and we have used example [T] with tti = 3/5,7r2 = 2/5 
and initial measure /tq in case a). 


To quantify convergence, we compute the difference |6j(/x„) — 6j(0„)| and we plot it in double 
logarithmic scale in fig. |9l We observe that the decay of this difference is bounded from above 
by An! 3 , in accordance with [5]. At large values of j this decay stops and difference between 
bjifJ-n) and hj{9n) oscillates at a low value due to the combination of two effects: on the one 
hand, computations are performed with a finite set of frequencies /C/,. When n increases the 
cardinality of /Cl increases as L", so that we are forced to reduce L in order to keep the number 
of frequencies within acceptable limits (see Appendix). On the other hand, our Fourier analysis 
is performed on finite Jacobi matrices whose length is fixed, albeit large (typically, 800,000). 
Therefore, amplitudes and phases of wJJ are determined less precisely when n grows. 

Next, let us consider case b). Figure fTOl is the analogue of fig. [9l but it now compares cases a) 
and b) at the fixed value of n = 2. The exponential decay of the latter is immediately appreciated. 
Plotting the same data on a larger interval, and resorting to logarithmic scale, shows the 0(l/j) 
bound on the decay: see fig. [TT] The exponential convergence of case b) appears even more 
dramatically here. The same comments on precision of the computation described above also 
apply in this case. 

We can now conclude this section noting that the data presented confirm the validity of our 
approach in determining numerically the isospectral torus for finite collections of intervals. We 
can therefore move to the Cantor set case. 


8 The isospectral torus of infinite gap sets—Cantor sets 

A recent paper by Kruger and Simon [39) discusses in the usual illuminating way various issues 
related to the definition(s) of an isospectral “torus” associated with measures supported on 
Cantor sets. Our approach to this problem is the following: in accordance with definition [T] we 
would like to define this object as 

Definition 2 The isospectral torus of a hyperbolic IPS Cantor set with finitely many maps is 
the set of Jacobi matrices associated with measures 9oo, defined as weak limits of 9n, when this 
latter is defined in eqs. U5116\ 17\) and A" is defined in eq. (0). 
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Figure 9: Difference between the outdiagonal entries of the Jacobi matrix of Hn and of its limit 
9n in the isospectral torus. Here, n = 2 (blue line), n = 3 (green) and n = 4 (red). We have used 
example [T] and in case a). 



Figure 10; Difference between the outdiagonal entries of the Jacobi matrix of /i„ and of its limit 
6 n in the isospectral torus. Here, n = 2 and fiQ is case a) (open red circles) and b) (filled blue 
circles). The curves plotted as reference (in reverse color coding to improve visibility) are an 
exponential (red) and a power-law decay A/j (blue). 
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Figure 11: Same data as in fig. [TOl now in double logarithmic scale. 


It is clear that the above definition can be extended beyond the realm of I.F.S. with a finite 
number of maps. Yet, we have been intentionally conservative, postponing generalizations. In 
fact, we believe that these hypotheses should be enough to assure weak convergence of the 
following procedure. Write the gaps in if" as: 

G" =G"-ilJi?", (22) 

where is the set of “new” gaps created at level n. is iteratively constructed as 

M 

i?" = U (23) 

i=i 

starting from = (/3i, a^) U ■ • ■ U(/^m- 1 ! c^m)- The countable set of gaps G°° in the I.F.S. 
Cantor set is then written as 

OO 

G°“ = y if". (24) 

n—1 

In so doing, we can order the countable set of gaps of the Cantor set as they appear with n in 
eq. (m, and for each iJ" in ascending order on the real line. This also imply an ordering of the 
infinite set of points {^i}ieN (one in each gap) and of the Riemann surface indices {(TijigN- 
Let us now consider such a double sequence that defines a sequence of mea¬ 

sures {0n}igN- If this sequence converges weakly, then the Jacobi matrix elements aj(0„) and 
bj{9n) also converge. As mentioned, we conjecture that IFS hyperbolicity should warrant this 
convergence. 

Our technique permits to test this conjecture numerically. As a matter of fact, convergence 
seems to hold at a geometrical rate—as it might be expected, because of the hyperbolicity 
hypothesis. In fig. [T2]we plot the absolute value of the difference bj(dn) — &j(0„_i), versus both 
n and j. On the surface defined by these values, we also plot the line N{e,n) indicating that 
for j < N{e,n) the difference between bjiOn) and its infinite limit is less than e. This quantity 
is significant, because it indicates that, at precision e, the truncated isospectral Jacobi matrix of 
rank N(e,n) is insensitive to the torus variables {Ci,cri}i>Q, with Q = M". As a consequence of 
geometric convergence, for any e, N(e, n) grows asymptotically in an exponential way with n. 
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Figure 12: Computing the isospectral torus of an infinite-gap set: see text for details. 



n = 2 

n = 3 

n = 4 

n = 5 

C? 

1.55434055 

2.4297566IE-002 

1.55445514 
2.44I07202E-002 

I.55438I55 

2.44958535E-002 

1.55429495 

2.45456487E-002 


Table 1: Fundamental frequency w” and related amplitude C" for the second gap in FI”, for the 
IFS of example [H 


The same geometric convergence is observed in the harmonic analysis of these Jacobi ma¬ 
trices. This holds for both frequencies and amplitudes. For example, in Table [T] we report the 
fundamental frequency corresponding to the second gap in our ordering, and its amplitude. This 
geometrical convergence might be the key by which to resolve the conjecture mentioned in the 
introduction, at least numerically. We now turn to this last point in this paper. 


9 If they are almost—periodic? 

0 As mentioned, the numerical techniques presented in this paper aim at providing hints to the 
resolution of the conjecture on the almost periodicity of hyperbolic IFS Jacobi matrices. We 
want to report here preliminary results in this direction. These results are not resolutive in a 
sense or another, but clarify the setting of the problem. 

For sake of definiteness, let Hoo be the balanced measure of such an IFS, its Jacobi matrix, 
and the sequence of its out-diagonal entries. Let also /i„ be the sequence of its IFS 

approximations, constructed starting from the Chebyshev measure of second kind, and the 
corresponding measures in the finite dimensional isospectral torus. 

We have seen that as j grows, bj{9n) and bj^Hn) approach each other exponentially fast. It is 
now important to gauge the slope of this exponential. In figure[T3]we plot the cases n = 1 to n = 4 
for the example considered in this paper. The fitting exponential lines have form They 

have been computed as a log-linear least squares fit, and appear to be approximately intersecting 
at n ~ I. The important feature is then the decay coefficient (i„. In table [2] we report these 
values. It appears that, as n grows, dn tends to zero. As a matter of fact, the sequence of 

^The title of this section was clearly derived from [ni¬ 


ls 











Figure 13; Differences \bj{On) — bj{ijLn)\ versus j, for n = 1,..., 4, and fitting exponential curves. 



n = 1 

n = 2 

n = 3 

n = 4 

dn 

0.300994 

0.164903 

0.080720 

0.044511 


Table 2: Coefficients of the fitting exponentials in fig. [13] 


dn values can be itself extremely well approximated by an exponential decay: dn e with 
S ~ 0.644. 

Let us now fix a positive threshold e. Extrapolating the data so far obtained, we can predict 
that for j > J{e,n) = log(l/£)e‘^"' the distance between bj{0n) and bj{ij.n) is less than e. 

Let us now consider the distance between and bj{^oo)- At difference with the previous, 

this distance can be bound by a fixed value e when j is smaller than a value N (e, n). We therefore 
define N{e,n) := max{/ s.t. — bj{fj,oo)\ < £,0 < j < 1}. It was originally found in [46] 

that N{e, n) increases exponentially with n. Notice that the same dependence has been observed 
in Sect. [8] and fig. [121 but we are now in a completely different setting, because one of the 
terms, fj,oc, is the I.F.S. singular continuous measure that is the object of the debated conjecture. 
Suppose that it happens that for e sufficiently small, N{e,n) ^ Then, if the inequality 
K> 5 holds, by choosing larger and larger values of n, the Jacobi matrix elements bj{^oo) can be 
uniformly approximated by bjlOn) in the asymptotically long, asymptotically remote segments 
j G (e'*",e«"). 

Our original technique [43) of computing the Jacobi matrix associated with /Too can be now 
put to work, in combination with that in [46] . to compute the function |6_,(/in) ~ &j(/ioo)|, that 
is plotted logarithmic scale in fig. [Td] At the basis, the experimental curves j{n) = N{e,n) for 
various values of e are plotted, as well as k{n) = with 6 given by the extrapolation of the 
data in table [2] Even if the figure spans over various orders n and up to Jacobi matrices of 
size in the hundreds of thousands, a clear-cut evidence in favor of the inequality k > 6 , or its 
violation, cannot be obtained: the displayed data seem to be very close to the hedge between the 
two behaviors. 
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Figure 14: Differences |6j(/r„) — bj{fioa)\ versus n and j. Also plotted at the basis of the picture 
are the values j{n) = N{s,n), for e = 10“^, 10“^, 10“^, 10“"^, and the (yellow) line k{n) = e"^", 
with 6 = 0.644. 

10 Appendix: Discrete Fourier Analysis 

The harmonic problem required in this paper differs from the usual in that the frequencies present 
in the signal are known, being determined by eq. (1201). Therefore, we act as follows. First, we 
multiply the “signal” bj by a window function Wj. For this, we choose the Dolph-Chebyshev 
window, well described in the classical paper |32] . Define the windowed discrete Fourier transform 

A (25) 

j 

Due to the property of the filter, the above is a finite summation. Next, we take into account the 
fact that the frequencies Wk are known. Since we restrict k to /Cl, this set is finite, and we can 
list frequencies in ascending order. Evaluating the windowed Fourier transform at any of these 
frequencies we find that 

F(c.k)= ^ lT(k,k')Dkq (26) 

k'eKL 

where Dk' is a complex coefficient that encodes both Ck' and '0k'- Actually, we work with two 
separate equations for real and imaginary part: this has several numerical advantages. 

kF(k, k') is a kernel that can be exactly computed, so that eq. (PSI) is a linear system for 
the unknowns Dk'. Taken at face value, this equation is discouraging, for the cardinality of the 
set K-l makes it untractable numerically. Yet, the Dolph-Chebyshev filter is such that VF is a 
function of the difference tuk — t-'-'k' > and this function can be made arbitrarily small outside of an 
interval centered at zero [32) . Therefore, in our numerical computations, by arranging frequencies 
in increasing order, fF(k, k') becomes a banded matrix, and the system (1^ can be swiftly solved. 
It is just the case to remark that this procedure yields results of a precision unattainable by the 
mere usage of the fast Fourier transform, that was applied to this problem in [44l |39] . 

Finally, a remark on the problem of extracting the phases ^k when examining the sequence 
bjlpn)- In this case, the finite size window Wj is shifted along the sequence, by a finite lag 1. A 
sequence of lags Im, m = 1,..., M is performed, and the corresponding amplitudes and phases 
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extracted. A fitting/extrapolation procedure is then performed that yields the asymptotic values 
Ck and V'k- The former are then compared with those obtained by the isospectral torus to assess 
the validity of the procedure, while the latter yield the desired, unknown phases. 
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